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REGULARIZED AND FRACTIONAL TAYLOR EXPANSIONS OF 
HOLDERIAN FUNCTIONS 

DIMITER PRODANOV 


Abstract. Holderian functions have strong non-linearities, which result in 
singularities in the derivatives. This manuscript presents several fractional- 
order Taylor expansions of Holderian functions around points of non- differ¬ 
entiability. These expansions are derived using the concept of a fractional 
velocity, which can be used to describe the singular behavior of derivatives 
and non-differentiable functions. Fractional velocity is defined as the limit of 
the difference quotient of the increment of a function and the difference of 
its argument raised to a fractional power. Fractional velocity can be used to 
regularize ordinary derivatives. To this end, it is possible to define regular¬ 
ized Taylor series and compound differential rules. In particular a compound 
differential rule for Holder 1/2 functions is demonstrated. The expansion is 
presented using the auxiliary concept of fractional co-variation of functions. 
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1. Introduction 

Difference quotients of functions of fractional order have been considered initially 
by du Bois-Reymond [9| and later by Faber uni in their studies of the point-wise 
differentiability of functions. The concept implies what is now known as Holder- 
continuity of the function. Later Cherbit [S] introduced the notion of fractional 
velocity as the limit of the fractional difference quotient. His main application was 
the study of fractal phenomena and physical processes for which the instantaneous 
velocity was not well defined. Independently, applications to physical systems ex¬ 
hibiting fractal behavior have inspired the development of local integral definitions 
of fractional derivatives. The definition of local fractional derivative introduced 
by Kolwankar and Gangal pT2] is based on the localization of Riemann-Liouville 
fractional derivatives towards a particular point of interest. Manipulations of this 
construct in many occasions are complicated and there are only few physical ap¬ 
plications so far. On the other hand, Ben Adda and Cresson mm demonstrated 
equivalence of so-defined local fractional derivative with an alternative definition 
based on a limit of a certain differential quotient. The correspondences between 
the integral and the quotient difference approaches have been further investigated 
by Chen et al. 5] and some of the results stated initially by Ben Adda and Cresson 
have been clarified recently [3] - 

The present work builds on the fractional variation operators, which in limit are 
equivalent to the fractional velocity introduced by Cherbit [6j . Previous work has 
demonstrated some basic algebraic properties of these operators and established 
their limiting behavior |16j . 
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The main application is a formal regularization procedure for a derivative of a 
Holderian function, which allows for removal of weak singularities in the derivative 
caused by strong non-linearities. This becomes possible by exploiting some of the 
properties of the fractional variation and velocity. One of main results of this 
procedure is the Taylor expansion of a Holderian function, termed the Taylor - ltd 
expansion, which makes use of a regularization of the integer-order derivative. This 
can be stated as follows: For a compound function f(w, x) where w(x) is H 1 / 2 the 
following expansion holds: 


d± s. df df d ± Id 2 f ± 

= S + S' W{X) ± 2S^ [W ' W] 


0 ) . 


The approach presented here can be used to characterize and approximate func¬ 
tions about points for which integer-order derivatives are not defined. Such func¬ 
tions are considered in the physical theory of scale relativity mm, mi where 
geodesic trajectories are considered to be of order H 1 / 2 . 

The manuscript is organized as follows: Section [5] gives general definitions and 
notational conventions. Section [3] introduces fractional (fractal) variation and frac¬ 
tional velocity. Section [4] demonstrates fractional and mixed-order Taylor expan¬ 
sions. Section [5] presents ordinary derivative regularization and demonstrates the 
compound differential rule for H 1 / 2 functions. 


2. Preliminaries 

The term function denotes a mapping / : R H- R or in some cases Ci—>C. The 
notation f(x) is used to refer to the value of the function at the point x. The term 
operator denotes the mapping from functional expressions to functional expressions. 
The symbol C° denotes the class of continuous functions. The symbol C n - the 
class of n-times differentiable functions where n £ N. Square brackets are used 
for the arguments of operators, while round brackets are used for the arguments of 
functions. Dom[f] denotes the domain of definition of the function f(x). 

Definition 1. Let i“ be the class of Holder continuous functions of degree a, 
a £ (0, i]. That is, we say that f is of class HI“ if V/ £ H a there exist two 
non-negative constants C, S £ K., such that for two values x,y £ Dom[f ], for which 
\x — y\ < 6, x ^ y, the following inequality holds 

\f(x)~f{y)\ <C\x-y\ a . 

For mixed orders n + a> 1, n£N, the following extended definition holds: Holder 
class HI n+a designates the class of C " functions for which the inequality 

I/O) - f{y) - Pn 0 -y)\<c\x- y\ n+a , 
holds in the interval [x, y\. P n {-) designates a real-valued polynomial of degree n of 

n 

the form P n (z) = &kZ k , where Pq(z) = 0 and a £ ( 0 , 1 ). 
fc=l 

Remark 1. Throughout the rest of the paper it will be assumed that, unless stated 
otherwise, the Holder exponent characterizing the function is point-wise unique. 
That is, there is only one fractional exponent associated with a point belonging to 
the domain of the function. 

In the opposite case we will discuss functions characterized by a Holder spec¬ 
trum, that is the set of such exponents. Functions having more than one exponent 
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will be denoted to as multi-Holderian and if necessary designated by the notation 
H“i + “2+- - . Indeed, it can be shown that if f(x) £ H Q then fix) £ H “ +/3 if /3 > a 
implying that H “ C H a+/3 . 

Therefore, unless stated explicitly the assignment fix) £ H Q will not denote 
sums of functions. 

Definition 2. Let the parametrized difference operators acting on a function fix) 
be defined in the following way 

A t [/] 0*0 : = fix + e) - fix) , 

A e “ [/] 0) := fix) - f{x - e) , 

A, [/] (x) := fix + e) - 2 fix) + fix - e) , 

where e > 0. The first one we refer to as forward difference operator, the second one 
we refer to as backward difference operator and the third one as 2 nd order difference 
operator. 


3. Fractional (fractal) variation and fractional velocity of 

FUNCTIONS 

Fractional (fractal) variation operators have been introduced in a previous work 
|T6l as: 


Definition 3. Let the Fractal Variation operators of order ft be defined as 


v [/] (*) := 

v e f[f] ix) := 


A+ If] (x) _ fjx + e) - fjx) 
e 0 e 0 
A £ - [/] (z) __ fix) - fix - e) 


where e > 0 and 0 < /3 < 1 are real parameters and fix) is a function. 


(1) 

( 2 ) 


In this paper we will consider the following extension of the previous definition 

Definition 4. Let the fractional variation operators of mixed order n + /3 acting 
on a function fix) be defined as 

f(x + e)-T n {x,e) 


<+/3 [/] 0 ) : = (»+!)! 

<+0 [/] 0) : = (-!)”(«+ !)! 


£ n+/3 

T n jx, —e) - fjx - e) 

,n+/3 


(3) 

(4) 


where T n (x, e), n £ N is the usual Taylor polynomial 

ta.ivixi+eA^c* 


fc=1 


k\ 


and e > 0 and 0 < ft < 1 are real parameters. 


Definition 5 (Mixed order velocity). The fractional velocity of mixed order n + f3 
of function fix) £ C n be defined as 

, n+ 0 f (nr\ -R™ , [J] (x) , 


off f ix) := lim v~ 

+ J V ' e—>0 


V- fix) := l™ v n+[) [/] i x ) 


(5) 

( 6 ) 


where e > 0 and 0 < a < 1 are real parameters. 
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Remark 2. The mixed-order definition obviously applies to functions that are in 
general differentiable, i.e. they may be non differentiable only at isolated points. 
In the case of non-differentiable functions, the fractional variation should be taken 
only in strictly fractional order. 

Theorem 1 (Modular variation). Let f(x) € C n . Then the fractional velocity of 
mixed order n + (3 about x can be evaluated as 


, n +P 


fi x ) = 


(n + 1)! Jim ^ ^ 


( n + P)(n- 1) £ ^° 


(n + 1)! Ar [/<")] (a) 


v n +h f ^ = v " ' ~ r — lim ■ 

(n +/3)( n _i) el 


(7) 


( 8 ) 


where e > 0 and 0 < /3 < 1 are real parameters and (n + /3)( n ) is the decreasing 
Pochhammer symbol, 

(n + P){k) = (n + /3)...(n-k + 0) . 

Let’s consider the point x € [x, x + e]. 

Proof. Forward case: The equation can be rearranged in the form 

/( x + e) - T n (x, e) = R n {x) , 

where R n {x) is the Taylor reminder. 

Then we have to evaluate the limit 


lim 


R n (x) 


e-S-0 e n+IS 

The evaluation of the limit can be performed by applying the l’Hopital’s 
rule n times. The Taylor reminder can be written in the form: 

x+e 

d n+l 


Rnix) = h J ( e + x ~ t)n fir+ 


Then 


n r x+e 

<«"<*». = a / 


d n+1 _i 

fjfTTlf (*) i x ~t + e) n dt = n R n - i(x) , 


where the index change in R n does not affect the derivative of fit). 
Therefore, the integral can be evaluated by induction to 


rx+e j n +1 


Ri = 


dt n+1 


f(t)dt=fW(t) 


x+e 


= At 


/ (n) l (*) 


The denominator can be evaluated by n-times differentiation about e in a 
similar way to give 

(e n+ 0!" ) = in + 0){n - 1 + /?)... (1 + = (n + %_ x) 

where we denote conveniently the decreasing Pochhammer symbol, x^ n ) ’■= 
x (x — 1)... {x — n + 1). Therefore, finally 
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Backward case: The equation can be rearranged in the form 
T n (x, -e) - f(x - e) = -R n (x) . 

The Taylor reminder can be written in the form: 

x—e 

d n+l 


R n(x) =-^ J {-e + x-t) r 


dt n+1 


f (t) dt 


Then we have to evaluate the limit 

i. R'U {x) 
inn- -5- . 

e ->0 e n+/3 

Applying again l’Hopital’s rule to the numerator yields 
[R n (x)]' e = -n Rn- i(x) 
using the same convention. Therefore, 

[Rnix)}^ = -(-l) n A 7 [/<">] (:r) . 

The expression for the denominator is identical to the forward case. 


□ 


Proposition 1. For a function f(x ) £ all mixed order fractional velocities 

exist up to the order n + p. 

Indeed, the proposition follows directly from Definition [T] for mixed order Holder 
functions. 


Corollary 1 (Continuous variation). If f( n \x) is continuous in the interval x £ 
[x,x + e] and f( n+1 \x) is continuous in the interval x £ (x,x + e] then 


v+ +fi f ( x) 
V- +0 f {x) 


(n + 1)! 

( n + /3)(n) 

(n + 1 )! 

( n + /3)(n) 


lim e 1 -^ f( n+1 \x + e) , 

e —^0 

lim e 1 -^ f( n+1 \x — e) . 


where (n + /?)(„) denotes the decreasing Pochhammer symbol. 


(9) 

( 10 ) 


Proof. Previous work has established |T6] that if the ordinary derivative f'(x ) is 
continuous in the interval x £ (x, x + e] then the following equations hold 


v+f(x) = i lime 1 0 f'(x + e) , (11) 

v-f (x) = ~ linr e 1-/3 /'(.x - e) . (12) 

p e—>0 

Application of PHopital’s rule to Theorem [T| demonstrates the desired result. □ 

Corollary 2. Let f{x) £ C” in the interval x £ [x, x + e]. Then the non-integer 
order velocities below the critical order n + P vanish. 
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Proof. Let’s calculate the fractional velocity of the mixed order m + 1 < n. By 
Corollary [T| 


m-\-/3 


f(x) = 


('m + 1 )! 


lim e 1 

e-K) 




0 m + a )(m ) 

where (3 < 1. Since by hypothesis f( n \x) is hnite everywhere in the interval the 
last limit is zero. Therefore all fractional velocities of order m + /3 < n vanish. The 
same argument can be applied to the backward fractional velocity. □ 


Remark 3. Following the notation of 0. 0 we define the usual forward and 
backward differential maps as: 

d + [/1W:= Z(£±fLTM 
Am w : = MzZlirf) 

These operators can be considered as particular cases of the Fractal variation ac¬ 
cording to Def. [3J 

Theorem 2 (Fractional Rolle). Let f[x) £ in [a, b] and f{a) = f(b) then there 
exists a number c £ [a, b\ such that v^ + [f] (c) < 0 and [/] (c) > 0 . Respectively, 
v^_f (c) < 0 and v+f (c) > 0 . If both variations agree then v®_f (c) = v+f (c) = 0 . 

Proof. The proof of the theorem follows closely the proof of the generalized Rollc’s 
Theorem. We can distinguish two cases: 

End-interval case: Since /(x) is continuous on [a, b], it follows from the 
continuity property that /(x) attains a maximum M at some ci £ [a, b] 
and a minimum m at some C 2 £ [a, b\. 

Suppose ci and C 2 are both endpoints of [a, b\. Because f{a) = f(b) it 
follows that m = M and so /(x) is constant on [a, b\. Hence [/] (x) = 
Vp~ [/] (x) = 0, Vx £ [a, b}. 

Interior case: Suppose then that the maximum is obtained at an interior 
point c of (a, b). For an e > 0 , such that c + e £ [a,b] f{c + e) < /(c) by 
assumption. Therefore, Vg + [/] (c) < 0. Therefore, v+f (c) < 0 . Similarly, 
for c — e £ [a, b] we have A 7 [/] (c) > 0 and v e p~ [f] (c) > 0. Therefore, 
v^_f (c) > 0 . Finally, if v^_f (c) = v+f (c) then both equal 0. The proof for 
the minimum follows identical reasoning. 

□ 


3.1. Bounds of variation. 

Lemma 1 (Bounds of forward variation). Let /(x) £ H“ in [x,x + e], v+f (z) ^ 0 

and C x and C' x be constants such that i) C is the smallest number for which 

|A+[/](x)| < still holds, that is C x = inf{|A/[/](x)| < Ceand ii) C' x 

c 

is the larqest number C for which |AjTfl(x)| > Ce@ still holds, that is C' = 

sup{|A+[/](x)| > Ce 13 } then 
c 

\vlf{x)\ = C x = C' x 

Remark 4. We will assume that v+f (x) is bounded in [x, x + e]. 
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Proof. We fix x so that C x remains constant. Let v+f ( x ) = a for a certain bounded 
a. From the definition of a Holder function it follows that |A+[/](:r)| < C x e^ then 
we subtract ae 13 from both sides of the inequality to obtain 

(C'x < |A+[/](ar)| - ae 0 < (C x - a) e 13 . 

Further division by the positive quantity e@ results in 


We have to further consider two cases. If A+[/](x) > 0 then taking the limit both 
sides gives 

Hm = „ 

e —>-0 

and lim C x — a > 0 but since C x is an infimum then lim C x = a. Since C is constant 
€—>-0 £—>-0 

with regard to e it follows that C x = a. On the other hand, if A+[/](x) < 0 then 

lim = _2a 
e->o eP 

and lim C x — a > —2a, therefore lim C x > —a and the same reasoning as in the 

previous case applies to yield C x = —a. Therefore, finally C x = |a|. The case for 
C' x can be derived using identical reasoning. □ 


Lemma 2 (Bounds of backward variation). Let f(x) gi 01 in [x, x+e], v ,3 _f (x) ^ 0 
and C x and C x be constants such that i) C is the smallest number C for which 
\Aj[f)(x)\ < Ce^ still holds, that is C x = inf{|A“[/](a;)| < Ce@} and ii) C' x 

is the larqest number C for which |Arffl(x)| > Ce@ still holds, that is C' = 

sup{|A“[/](x)| > Ce@} then 
C 

| vif(x)\ = C x = C' x 

Proposition 2. Let f{x) £ H“ about x and 0 < j3 < a < 1 then v+f (x) = 

v-f(x) = 0 . 

Lemma 3. Let f{x) £ H.P in the interval [a, b\. Then either v+f (a) ^ 0 or 
v+f (b) ^ 0 or v^_f (a) ^ 0 or v^_f ( b ) ^ 0. 

Proof. The proof follows from the Lebesgue’s differentiation TheoreuQ- The Holder 
condition fulfills the assumption of the theorem. Therefore, the condition of Corol¬ 
lary m eq. ED holds. Therefore, if /'(:r) diverges at a then v+f (a) ^ 0 holds. 
Similar considerations are applicable in the other cases. □ 


Theorem 3 (Discontinuous variation). Let f{x) £ H /3 in the interval [x,x + e). 
Then if v+f (x) ^ 0 then v+f (x) is discontinuous about x. If v^_f (x) ^ 0 then 
V-f {x) is discontinuous about x. 

Proof. Suppose that v+f (x) £ 1“ and v+f (x) ^ 0. 


1 Let f(x) : [a, b] —>• R be an monotone function. Then fix) is differentiable almost everywhere 
in [a, b]. For an elementary proof see m 
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Constant case: Suppose that v+f ix) = K ^ 0, constant. Then 


lim 

e-S-0 


f{x + e) - f(x) 
e? 


= K 


without loss of generality let’s assume that f(x a ) 
Therefore, 


lim 

e ->0 


f{x s + e) 
e? 


= K 


0 for a given fixed x s . 


and f(x + e) is of the form g(x s )e /3 . Therefore, g{x) is constant. Therefore, 
by the reasoning of the previous case we conclude that for /3 < 1 the 
hypothesis is false. 

Non-constant case: Suppose further that v+f (x) ^ 0. Then the conditions 
of Lemma [3] hold. According to Corollary []] the co-domain of v+f ( x ) 
or respectively v^_f [x) is discrete. Therefore, so-derived functions are 
discontinuous. 


□ 


4. Taylor expansions 

4.1. Fractional order Taylor expansions. 

Theorem 4 (Fractional Taylor expansion). Let fix) £ C° about x in the interval 
x £ [x,x + e]. Then the existence of v+f (x) at x for a (3 < 1 implies that 

f{x + e) = fix) + Ke p + 7 e ' 3 , 

where K = v+f ix) and lim 7=0 and vice versa. The existence ofv^_fix) implies 
that 

fix -e) = fix) - K*eP + 7 e p , 

where K* = v^_f (a:) and lim 7=0 and vice versa. 

e—>o 

Proof. Suppose that 

fix + e) = fix) + Ke p + 76 ^, 

where lim 7 = 0 . Then v e ^~ [/] ix) = K + 7. Taking the limit provides the result. 
Suppose that 

fix-e) = fix)-K*eP+ 7 ^, 

where lim 7 = 0 . Then —Vp~ [f] ix) = — AT* + 7. Taking the limit provides the result. 

□ 

Therefore, if follows that the following fractional Taylor expansion holds: 

fix ±e) = fix) ± vlf ix) e“ + O (e“+) , (13) 

where the last expression in the notation is interpreted as an arbitrary exponent 
larger than a. 
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4.2. Higher order fractional Taylor expansion. 

Definition 6. Let us define formal fractional Taylor polynomials for an increasing 
sequence {oi}" =1 as: 

n 

T a,n( X i e ) = f( x ) + 5Z C fc £ “ +fc > 

k—0 

n 

T a,n( x > e ) = f( x ) + 5Z c fc(-e) fe e Q+fc . 
fc =0 

where a denotes the multi-index and Ci are arbitrary constants. 

We will establish the relationship to the growth of the function at x ± e, respec¬ 
tively. We will look for O ( e m ) equivalence. That is 


f(x + e) -T+ (x,e) 

lim- : - = 0 

e-s-0 e m 

Obviously for k = 0 holds Co = vff(x). Then the second coefficient can be 
calculated as 

f(x + e) - f{x) - v%f (x) e a T+ (x, e) - f(x) - v%f (x) e a 

lim-—-= lim —- 

€—^0 e m e —>-0 6 m 


= lim 

e-X) 


Cie a+1 + 0(e 2+a ) 


e" 

= lim Cl e a+1 - m + lim 0(e 2+a ~ m ) 

e-s-0 e—SO 

Therefore, in order for the RHS to be finite we must have m = 1 + a. Then for the 
LHS we have 


lim 

e-s-0 


f( x + e) - f{x) - y+f ( X )! 

f l+a 


1 


lim 

1 + a e-s-0 


f'{x + e)-a vlf ( x) I 


The limit can be evaluated by application of L’Hopital’s rule and rationalization. 


1 r f'(x + e)-avf_f{x)e 

lim 


at —1 


1 + a e-¥0 


1 ef'(x + e)-avf_f(x)e° 

lim ■ 


1 + a e-sO 
lim ■ 


d+ot 


X f{x + e) - a v%f(x) 


1 + a e-S-0 


lim S- e 1 a f'(x + e) 


Therefore, 


1 + a e^o de 

ci = -p—— lim e 1 ~ a f'(x + e) . 


1 + a e-so de 

The same procedure can be extended for the general case by induction. For an 
arbitrary k < n we will have 

.. /(z + e) — T+ (x,e) f(x + e)-T+ k (x,e) 

Cb = lim-=—-= hm ■ 

e-s-0 e k+a e— s-o 

Therefore, we would have 


pfe+a 


1 fi x + c) - (r£ k (x,ej) 


lim 
k + a >o 


ck+a 
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1 


i f'(x + e)—^2 c j(a + j)e a+j 1 1 “ ^ 


lim ■ 

k + a e—>o 


i =o 


fk+a 


1 


,1 — O' i*/ 


lim ■ 


fix + e)-J2c k {a + k)e k 
1 


By applying k times L’Hopital’s rule the denominator can be evaluated to give fc! 
in order to eliminate the Taylor polynomial. Therefore, finally 

1 F) k 

- 6 e 1 -«/'(*+ e). 


Ck — 


lim 


/c! (1 + a) e—<9e fe 
Applying similar procedure to the backward case would yield 


Theorem 5 (Mixed order Taylor expansion). Let fix) € C n about x in the interval 
xG (x,x + e] or [a; — e,x), respectively. Moreover, let v+f (x) or vff(x) exist 
and are non zero at x. Then the following expansions hold: 


fix + e) = f(x) + c ^ a+k + O (e n+a ) , 
fc—o 


where 


and 


where 


c k = TTTi -r lim T—r e 1 a f'(x + e) 

k\ ( k + a) e->o de k v ' 

n 

fix - e) = fix) + J2 Ck * a+k + ° ( £n+Q ) 

fc =0 


Cfc 


(_l)fc gk 

k\ ( k + a) «->o de k 


1—a 


fix 


e). 


Example 1. Let’s compute the fractional velocity of f(x) 
x = 0. The first derivative of the function is 


fix) 


1 

\/2 x — x 2 ’ 


arcsin(l — x) about 


which is undefined about 0. Therefore, the function does not posses a Taylor ex¬ 
pansion about 0. On the other hand, 


v+f (°) 


— — lim 


13 e-s-o v/2 


2 


— e 


Therefore, for (3 = 1/2 v+f ( 0) = — yfl. The fractional Taylor expansion about 0 
then is 

fi x ) = ^ - v /2 y/x + O(fx) 

The regularized derivative then is 


df(0 ) = lim ■ 

e ->0 


arcsin(l 


e) — — + y/2 y/e 
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Figure 1. Error of the mixed-order Taylor expansion of arcsin(l — 
x) about the origin 


lim Z-(Z-- —L= 

e_> ' 0 v^Vv^ a /2 - e 


= — lim ■ 




e_K) a/2 - e 

Applying Theorem^ gives the following approximation: 


= 0 


f(x) « x - V2^x~ — 3 - 

^ O ■ 2 2 


O ^ 

3 • X 2 


5 • X2 


_ -.a _ ^ 

5-22 7 • 2 T 


35 • a; 2 
9 . 


The approximation can be appreciated in Fig. [7] where the error of this expansion 
is plotted. 


5. Derivative regularization 

5.1. Fractional order regularization. Here we will repeat without proof results 
that will be used further in the regularization procedure: 

Proposition 3. Let /( x) £ C 1 about x and 0 < /? < 1 then v+f (a;) = vH_f (x) = 0. 
Proposition 4 (Duality of the limit of variation). If /( x) is defined in [x — e, x + e], 

the limit lim [/] ( x ) exists and lim — 6 ^ ^ ^ = 0 then v+f ( x ) = v^_f (x). The 

proof is detailed in |16] . 


The value of the fractional velocity of a function f(x) of fractional grade about 
a point can be regularized to zero using the following result: 


Theorem 6 (Fractional variation regularization). Let f(x) £ C° and v+f (x) 
and/or v 0 _f (x) be defined at x for a given 0 < /3 < 1 then 


where g(x, e) 
kernels. 


lim f{X + e) 7 9{X ' e) = lim »•(*■<>-/(*-«) 

e -)-0 6 ^ e —>-0 

f(x) + v 0 + f( x) e 0 and g*{x, e) = f(x) — v 0 _f (x) 


= 0 

e 0 are regularizing 
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Proof. Let g(x , e) be smooth about the argument e. We will consider further four 
cases. 

Differentiable /( x): Forward case : Let /(x) £ C 1 in the left-open interval 
(x,x + e] and g(x : y) £ C 1 by its second variable in the closed interval 
y £ [0, e]. Then we can apply l’Hopital’s rule to the limit 


lim 

£->0 


f(x + e)-g(x,e) 


\iraj——{f'{x + e)-g' e {x,e)) 

e->0 p 

v+f( x ) - \ lim e 1_/5 ff'(ar, e) 

P £->0 


by Th. |T] Let’s suppose that je 1 ^'(x, e) = K. Then the resulting 
ordinary differential equation can be solved to yield for g{x, e): 

g(x, e) = C + Ke& . 

Therefore, K = v+f (x). Moreover, setting C = f{x) allows for interpola¬ 
tion with the fractional variation operator. 

Backward case: 

g(x, e) - f[x - e) e 1_/3 , ,,, u \\ 

h n 0 - j -= iTo ~t (/ ( _ e)+5e( ’ e)) 

= V p _f{x) + J lim e 1 ~ 0 g' e (x,e) 

P £-10 

by Th. [i] Therefore, in this case K = —v^_f (x). 

Non-differentiable f(x): Forward case: According to Th. [4] the existence 
of v+f (x) implies that 

f(x + e) = f(x) + Ke p + 7 , 

where K = v+f (x) and lim 7 = 0 . Therefore, 


f(x + e)-g(x,e) f{x) + Ke p + pe p ~g{x,e) 

lim - o -= lim —— - o - 

e—>-0 e—>0 


= lim 

e—>0 


f{x) +KeP -g(x,e) 


Therefore, since the value f{x) is arbitrary we can set g{x, e) = f{x) + Ke^. 
Backward case: 


f{x - e) = fix) - Ke@ + 7 
where K = v+f (x) and lim 7 = 0 . Therefore, 


£—10 


lim gjx,e) - fjx-e) = lim - fix) + Ke p - 7 + g(x, e ) 


£-10 


£— 10 e@ 

-fix) + Ke 13 + g{x,e) 


= lim 
£-10 


Therefore, in this case gix, e) = /(x) — Ke@. 


□ 
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Remark 5. The proof can be extended also for /3 = 1 (Differentiable case) where 
it can be noted that <?'( x,e) = K and therefore g(x,e) = K e + C implying that 
9 (x,e) = f(x) + f'(x) e. 

The last Theorem can be also used to extend the usual l’Hopital - Bernoulli rule 
for the indefinite form 0 / 0 : 

Theorem 7 (Fractional l’Hopital Bernoulli, case 0/0). Let f(x) = 0 and g(x) = 0 
for a given x for two functions /, g £ C° and v+g ( x ) ^ 0 identically in [x, x + e], 
Then 

v+f{x) r f(x + e) 

-XL -= lim —-- 

v+g (x) e_>0 g(x + e) 

v-f(x) r f(x-e) 

v_g(x) £_>0 g{x — e) 

Proof. Forward case: By assumption f(x + e) = f(x) + Ke@ + ye * 3 = Ke@ + 
setting v+f (x) = K. g(x + e ) = g(x) + Pe ^= Pe^+ 7 *e^ setting v+g ( x) = P. 
Then we have 

f(x + e) _ K + 7 
g(x + e) P + 7 * 

Taking the limit gives 

v f{x + e) K vlf(x) 
e ^° g( x + e) P v+g (x) 

The same reasoning can be applied for the backward case. □ 


Definition 7 (Composition of variations). Let f(x) be continuous about the point 
x. Let the composition of fractal variations be defined as 

v A+[/](x)-u^/(x)e /3 

lim-;—- 

€->0 


a°/3/ ._ Jj m 


u“°' 3 / (x) := lim 

e-J-O 


e a+/3 

a e ~[/](x) -i//(x) e 13 

f a +p 


(14) 

(15) 


where 0 < a, /3 < 1. 

Remark 6. We will assume that v+f (x) = 0. Indeed, since f(x) is continuous 
about x then lim A+ [/](x)] = lim A“[/](x)] = 0. 

e —>-0 e —>-0 

Proposition 5. Under this definition for the general case a ^ /? it follows that 
v% o0 f{x)^vl oa f{x). 

Proof. Indeed if we take the difference 

A+[/](x) - v+fjx) e p _ A+[/](x) -v%f{x) e a _ v^f (x) e? - vff (x) e a 

e ct+,8 £ a+/3 g«+|S 

If we take the limit and apply L’Hopital’s rule for e we get 

1 


— lim 
e-K> a +, 


1—(a+/3) 


(M/( 


x) e' 3 1 — av[ 


-f(x)e a - 1 ) 


a + 

Therefore, if a ^ /3 the limit is diverging. 




□ 
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Definition 8 . Fractional order regularization is defined in the following way: 

d 0+ t , ^ 1—00/3 r A+[/](x)-^/(x)e^ 

——/(x) := u + M 7 (x) = lim--- 

ax e->o e 

d P ~ t , y 1 —/3o/3 j: ( \ y A 7[fK X )~ V -f{ X ) eP 

—f( x ):= v _P f (x) = lim--- 

We will require as usual that the forward and backward regularized derivatives be 
equal for a uniformly continuous function. We will use further also a shortened 
notation df(x) when the value of (3 is fixed. 


By direct calculation we can use the composition of variations to regularize the 
first derivative about non-differentiable points. The proof follows from Th. [ 6 ] 
Having established the bounds of variation we can proceed to a new result for 
a Holderian function of grade a. Under this regularization we have the following 
result: 


Theorem 8 (Orthogonality of grades). Let f £ HI", which is not a sum o/C 1 and 
H Q functions. Then 

d a f(x) v$f (x) = d a f(x) v°ff (x) = 0 , 

where we interpret regularization with regard to the grade a. 

Proof. Let /(x) £ H“. We fix x so that C x remains constant. 

C' x e a < |A+/(x)| < C x e a 

Let v+f ( x ) = a assuming that a is bounded. Subtracting \a\e a from both sides 
gives 

& - M)e“ < |A+/(x)| - |a|e“ < (C x - |a|)e“ 

Therefore, if we take the supremum and the inGmuni bounds on C' and C respec¬ 
tively we will get 

\ A tf( x )\ - W a = 0 

There are two cases to consider. If a > 0 implying also A+/(x) > 0 then 


A+/(x)-u“/(x)e“ = 0 
If a < 0 implying also A+/(x) < 0 then 

-A +/(*)+«“/(x) e“ = 0 

Therefore in both cases 

A+f(x) — r°/ (x) e a = Q 
e 

but this is the definition of the regularized derivative. Therefore, if follows that 

n c 


dx 


■f{x) = 0 . 


□ 


Example 2. A key aspect to the former proof was the requirement ot keep C x 
remains constant. We will see that its relaxation invalidates part of the result. Let 
fix) = x + y/x. Then v+ 2 f (0) = 1. 

A + [/] (e) — \[e = e 
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Therefore, 3/(0) = 1. Inspection of the equations shows that C x is a function of 
the increment e. Therefore, the hypothesis can not be relaxed. 

5.2. Mixed order regularization. The regularization procedure can be extended 
also towards mixed orders in the following way: 


Definition 9. Let us define formal fractional Taylor polynomials for an increasing 
positive sequence {oz}" =1 bounded from above by 1 as: 

Ti,(x, e ) = /(x) + J> 

Ot-i k= 1 


T a ,n( X ^) = f( X ) - E 


c*e ai + 


E 

k=1 


/ (fe) ( a 

k\ 


-(-*)* 


where a denotes the multi-index and Ci are arbitrary constants. Then a-regularized 
derivatives are defined as 


d n+ 


dxn -f{x) = {n+l)\\im 


f{x + e) -T atU (x,e) 


e n+P 

T a,n( X ^) ~ f( X ~ e ) 

f n+P 


where f3 = supo^. We will require as usually that the forward and backward regu¬ 
larized derivatives be equal for a uniformly continuous function. 


Proposition 6. For a multi-Holder function H“ 1+a2+ ' according to Th. O we 
have Ci = vff f ( x ) and c* = v°f f (x). 

The proof follows by induction considering that the sequence {a,} is increasing. 
Then according to Th. [ 6 ] a new regularization term is added for every Oj. 


Proposition 7. Let us consider the following fractional Taylor polynomials for a 
positive increasing regular sequence {oi }" =1 : 

n 

T a,n( x > e ) = f( x ) + E Cfee “ fe > 
k= 1 
n 

T- n (x,e) = f(x) + J2c k (-l) k e ak , 

fc=l 

where a denotes the multi-index and Ci are arbitrary constants. Then the following 
expansion holds for a compound function /( x a ), a 6 (0,1): 

n 

/(* + <0 = /(*)+£c*e a * + O( C n “) , 

fc=l 


for 


1 d h 

c * = t\ a# fiu) 


u—x a 


f( x -e) = f(x) + Y / c k e ak + 0(e na ) , 
k =1 


and 
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for 


_ (- 1 )* d k 
Ck k\ du k ^ 


Proof. We will establish the relationship to the growth of the function at x ± e, 
respectively. The technique of the proof is similar to the previous case. 

Obviously for k = 1 holds ci = v+f (x). Then the second coefficient can be 
calculated as 


lim 

£->0 


f(x + e) - f(x) - vlf (x) ( 


= lim 

e^O 

= lim 

e - K ) 


T a,n(x, e) - f{x) - V%f (x) i 

e m 
2 q\ 


cie 2a + 0(e 2a ) 


= lim cie' 

e—>0 


2a—m 


lim 0(e 

e —>-0 


2 a—m 


Therefore, in order for the RHS to be finite we must have m = 2a. Then for the 
LHS we have 


lim 

e-s -0 


fix + e) - f{x) - v^f jx) , 

/2 a 


1 f'ix + e)~av‘ffix)> 


= -lim 

2 a e— 


c 2a-l 


The argument of last limit is then 


c a—l 


, 1 —a 


c 2 a—1 


-ie^fix + e)^v^fix)) = 


f'ix + e)-vlfix) 


The limit can be evaluated by application of L’Hopital’s rule and rationalization. 


A—a fi 


lim • 

e->-0 


f'i x + e ) - v +f (x) 


= — lim e 1 
a e->o 


A 

de 


ie'-y'ix + e)) 


Therefore, 


c 2 = TT -2 !™ el a 4- 

2 or e^o de 


_l—a 


d_ 

We 


fi x + e) 


However, this can be recognized as double application of the ordinary compound 
differential rule. Indeed, let us suppose that 


fix) = fiuix)) , uix) = x° 


Then 


and 


duixf \ 1 dfjx) = dfjx) 
dx ) dx du 


9 fix) x 1 a dfix) 


du a dx 

and we can recognize the derivative of x a evaluated at x = e. Similar arguments 
can be demonstrated to hold for all k < n by induction. 

Therefore, in the general case we would have as expected expansion of the func¬ 
tion fix) in Taylor series w.r.t u(x) = x a followed by substitution. Therefore, 


for 


/(® + e) = /(*)+£c*e afc +0(e na ) 

/(«) 


fc=l 

1 d k 

k\ du k 


C-k — 
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and 


for 


f(x-e) = f(x) + J2cke ak + 0(e na ) 

k =1 




(—l) fc d k 
k\ du k 


f(u) ■ 


□ 


Remark 7. In the original notation we can also state the result explicitly 

\k k 


Ck = 1 ™!]' 


k\ 


) /(*±e) 


(16) 


where the product has to be understood as composition. 


Example 3. Consider the function f(x) = cos(t 1 / 3 ). We will develop the frac¬ 
tional Taylor expansion about x = 0. The first derivative of the function is 


/'(*) 



which is undefined for x = 0. According to Th. [5|/or a = ^ we would have 


Cl 


Continuing for other degrees 


— lim sine 1 / 3 = 0 

e-S-0 


Co — — ~ lim cos e 
2 £—>-0 



and so on. Therefore, in agreement with the previous proposition we would have 


f(x) 


, 1 2 1 

1 -• x 3 H- 

2 24 


4 

13 


1,1s 

-- x -\ -• X 3 + . . . 

720 40320 


6. Compound differential rules 

Theorem 9. Let Let f(t,w ) € C 1 be composition with w(x ) € H ^ in the interval 
[x, x + e] then the following expansions hold: 


p f ( \ df p , x 

V+f (W) = gff V + W ( X) 

p r ( X df p 

V_J (W) = -g — vLw ( x) 

Proof. According to m we have 

v P + [f}M = ^f^l v e+ [w]{x) . 

Taking the limit provides the result. The backward case is proven in a similar 
way. □ 
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6.1. Regularized Taylor — Ito expansion. The differentiation rule for compo¬ 
sitions of functions can be derived using the regularization of the derivatives. It 
provides a result formally analogous to ltd’s Lemma. In a similar way the rule can 
be derived using a non-differentiable embedding principle [ 8 j. 

Theorem 10 (Taylor - Ito expansion). Let Let f(t,w ) £ C 2 be composition with 
w{x) £ HI in the interval x £ [x, x + e\ then the following expansions hold: 


d + ,, , df df d+ Id 2 f 

si + ST nr w{x) + 2 dwt [w ' w] 


+ V) 


d df df 

S- f{t ’ w) =M + -^ 


d II 1 ° 2 f r 1 —/ \ 

dr"’ w -2s^ K “' 1 (l) 


Proof. We will prove first the forward case. 

Forward case: By second order Taylor’s expansion we get 


df 

f{w(x + e), x + e) =f(w, t) + e 


—A + 

dw e 


0*0 + 


i5V 2 

2 9x 2 6 


5 2 / 

dx dw 


(x)A+ [w] (a:) e+ 


15V 

2 9w 2 


(A+ [w] (V ) 2 +C>(e 2 ) 


It follows that two cases have to be considered: 

Differentiable case 

This case is realized when vV 2 w ( x ) = 0 and the ordinary derivative is 
defined: 


d t [/] K V =lr + 7P7 d t M 0*0 ■ 


dx dw 

15V 5 2 / 

2 9a; 2 6 5a; 9w 


19V 

2 9w 2 




(x)A+ [w] (x) + 

O(e) 


Taking the limit gives the expected result 

d£ _ 9/ 9/dw 
dx dx dw dx 


Non-Differentiable case 

1 /2 

This case is realized when n + 7 w (x) = K ^ 0 and the ordinary derivative 
is not defined: In this case w(x + e) = w(x) + (K + where lim 7 = 0. 

e—>-0 

We substitute partially in the Taylor expansion: 


/ (w{x + e), x + e) — / (w, x) = ^ £+^-(^ + 7 ) 2/6 

ax au> 


19V 

2 dx 2 

1 9 2 / 


£2+ sV l)A - +W(l)e+ 

(a+[»](x)) 2 + o(e 2 ) 


2 9u> 2 
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Then after rearrangement we get 


/ (w(x + e), x + e) - f (w, x) - ^-(K + 7 )^ 


df , ld 2 f 


^ ^ -(x)Af [w] (x) + - 77 -C [w,w;l + 


0(e). 


dx 2 dx 2 dx dw v c 2 dw 2 

But the LHS is the forward regularized derivative of /(w, x). Therefore, we 
obtain 


u+f _ df 1 d 2 f + 

dx dx + 2 dw 2 ^ ^ 

Therefore, we can rewrite the expansion using the uniform notation 
d + f df df dw 1 d 2 f , 

~dhT = dx + d^ &T + 2 d^ 2 
Backward case: By second order Taylor’s expansion we get 

f(w{x -e),x-e) = f{w , x)-^j-e- [w] (x) + 


(17) 

(18) 


1 

2dx 2 

1 d 2 f 


d 2 f 

dx dw 


(x)A e [w] (x) e+ 


2 dw 2 

Differentiable case 

1 /2 

When v_ w (x) = 0 and the ordinary derivative is defined: 

dr [/] (w, x) = ^ + |^d e - [in] (x) - 

d f (x) A 7 M (x) - iu KHr + °(t) ■ 


10V_ 

2 dx 2 6 dx dw 


2 dw 2 

Taking the limit gives the expected result 
df_ _ df_ dj^dw 

dx dx dw dx 

Non-Differentiable case 

when v_w ( x ) = K 0 and the ordinary derivative is not defined: 
In this case w[x — e) = w(x) — (Krfy/e where limy = 0. We substitute 

e—>-0 

partially in the Taylor expansion: 

/ (w, x) - / (w{x - e), x - e) = e - ^(K + 7 )Ve 

1 d 2 f 2 d2 f f A - r 1 ^ 

“2 W e d^ {x)A * [w]{x)t - 
1 d 2 f 


2dw 2 

Then after rearrangement we get 


(A+[w] (a;)) +0(e 2 ). 


f (w, x) - f (w(x -e),x- e) + §£{K + 7 )^ 


d£ 1 d 2 f 

dx 


2 dx 2 6 dx dw 


(x)A+ [w] (x) - [umo1 + 


2 dw 2 


0(e). 
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But the LHS is the backwards regularized derivative of f(w,x). Therefore, 
we obtain 

d~f df 1 d 2 f 

IT = Si ~ lw ’ w] ~ (19) 

Therefore, we can rewrite the expansion using the uniform notation 
d~ f df df dw 1 d 2 f 


dx dx dw dx 2 dw 2 


\w, w\ 


( 20 ) 

□ 


From the proof of Th. [TOland the results of Cresson and Greff [5] it follows that: 

Proposition 8. For a composite function f(t,w) £ C 2 where w(t) £ H 1 / 2 the 
following approximation holds 

df r df 1/2 ^ , e ^ d 2 f 


At [/] (x,w) = e — + V^g^v 1 + /2 w(x) + - J2 


2 dxj dxj 


[Xi, Xj]t + 0 (y/e) 


( 21 ) 


A e [/] (x, w ) = e ^-^^- vt /2 w nlLj Xh X ^ + °W ’ 


2 dxj dxj 


( 22 ) 


where e = dx. 


The last result can be used for numerical approximation of fractional singular 
problems. 

Example 4. Let’s calculate the expansion of a quadratic function f(v) = —v 2 {x) 
composing with a H 1 / 2 function v(x). 

^f(x,v) = v(x)^- + ^[v,v] + (x) 

For example, if v(x) = y/x we have for x ^ 0 

A /(l , = ^ _T + l lv ^, = I + ° = I, 

which is consistent with the derivative of f{x) = — x. For x = 0 we have 

A/(z,.) = V5-0+1 = 1, 

which is again consistent with usual analytical treatment. 


7. Discussion 

Classical physic variables, such as velocity or acceleration, are considered to 
be differentiable functions of position. On the other hand, quantum mechanical 
paths m. 0, mi and Brownian motion trajectories were found to be typically 
non-differentiablc. The relaxation of the differentiability assumption opens new 
avenues in describing physical phenomena, for example, using the scale relativity 
theory. 

The aim of the present work is partly to explore the mathematical foundations of 
the theory of scale relativity championed by Nottale E. H, E' The main tenet 
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of the theory is that there is no preferred scale of description for physical reality. 
Therefore, a physical phenomenon should be described simultaneously at many 
irreducible scales. This assumption naturally leads to fractal and non-differentiablc 
representation of the studied phenomena. 

Specifically, the derivative regularization procedure demonstrated here and the 
resulting Taylor - Ito expansion give support to some of the derivation procedures 
employed in scale relativity. 
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